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ABSTRACT 



An improved model for a once-through counter-cross-flow 
waste heat recovery unit with a segmented fin-tube arrange- 
ment was developed. The model was coded in FORTRAN IV 
computer language for use on an IBM-3033 computer system, 
and was tested for various conditions of uniform and non- 
uniform gas glow distributions. Additional parameters 
varied were heat exchanger size, steam outlet pressure, 
steam outlet temperature and gas flow rate. The effect of 
each variation on WHRU performance was evaluated and the 
interrelationship discussed. Results indicate the ideal 
design is a high steam temperature, low operating pressure, 
flow distribution controlled heat exchanger. 
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INTRODUCTION 



I . 



A. BACKGROUND 

Since the introduction of the DD-963 class ships into 
the fleet, the U.S. Navy has become firmly committed to gas 
turbine propulsion plants. Currently, besides the DD-963 
class, there are four additional ship classes which are, or 
will be, gas turbine powered. These are the Patrol Hydro- 
foil Missile (PHM's), the Oliver Hazard Perry FFG-7 class, 
the Ticonderoga CG-47 class and the soon to be built DDG-51 
class . 

This commitment to gas turbine propulsion has stimulated 
great interest in the concept of converting the waste heat 
from the gas turbine engines into useful shaft power. In 
December 1979 the Naval Sea Systems Command awarded contracts 
for development of a conceptual design of a new generation 
of waste heat recovery systems. The goal of the design was 
to provide a system which would realize a 25 percent minimum 
improvement in ship propulsion fuel consumption and increase 
the ship's operating range by 1000 nautical miles. The out- 
growth of the effort has been the RAnkine Cycle Energy Recovery 
(RACER) system. As conceived, the RACER system will be an 
unfired waste heat recovery system designed to convert waste 
heat from the exhaust of the main propulsion gas turbine 
engines into useful shaft horsepower. Such a system has 
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commonly been called a COGAS or Combined Gas And Steam system. 
The hot exhaust gases of the main propulsion turbine are 
passed through a Waste Heat Recovery Unit (WHRU) , where steam 
of the desired temperature and pressure is generated. The 
steam is used to drive a steam turbine which operates in 
parallel with the propulsion gas turbine through a common 
reduction gear. 

Numerous papers have been written concerning the RACER 
system. Halkola, Cambell and Jung [Ref. 1] discussed the 
conceptual design for the system. Mattson [Ref. 2] presented 
a plan for designing reliability and maintainability into 
the RACER system and Combs [Ref. 3] developed a computer model 
of a COGAS system which designs a waste heat recovery unit 
and allows comparison of system performance by variation of 
system parameters. 

References 1, 2 and 3 all agree that the desired goal of 
fuel conservation must not come as a result of increased 
maintenance, manpower, or main propulsion system down-time. 

The design philosophy dictates that the system have high 
reliability, be inherently safe, have low maintenance require- 
ments and be essentially automatic in operation. 

Reference 1 states that the waste heat recovery unit 
paces the design of the system. The majority of effort by 
Combs [Ref. 3] was to create a WHRU design model which would 
simulate the performance characteristics of a physical boiler. 
Since great expense is involved in the design, fabrication 
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and operation of prototype systems, it is of primary impor- 
tance to achieve a realistic mathematical simulation first. 

To this end. Combs devised a model which designs a WHRU 
when given feedwater inlet temperature, steam outlet tempera- 
ture, operating pressure, inlet and outlet gas conditions 
and gas flow rate. 

In this study, the initial model of Combs was considered 
and an improved computer model of the waste heat recovery 
unit was developed. A once-through counter-cross-flow heat 
exchanger was selected because of the advantages of lower 
initial cost, faster response, compactness and low operating 
cost. The improvement in the model is that the heat exchanger 
is discretized into individual "nodes" much in the manner 
done by Shu and Kuo [Ref. 4], This procedure treats smaller 
physical pieces of the overall heat exchanger as individual 
small heat exchangers to attain a better understanding of 
the temperature distributions within the WHRU. The improved 
model allows an analysis of the effect of non-uniform velocity 
distributions on the performance of the WHRU, and ultimately, 
the RACER system. 

B. OBJECTIVES 

There were five main objectives for this thesis. 

1. Devise an improved computer model for the once-through 
WHRU. This model was to predict the performance of various 
size WHRU's given feedwater inlet temperature, steam outlet 
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temperature, WHRU operating pressure, inlet gas mass flow 
distribution and inlet gas temperature distribution. 

2. Enable the designer to establish a design WHRU size 
based on the performance of the model with various input 
conditions . 

3. Evaluate the design WHRU utilizing various uniform 
and non-uniform inlet gas distributions for both on-design 
and off-design conditions, 

4. Re-evaluate the WHRU utilizing the same uniform and 
non-uniform gas distributions into the WHRU, but with 
different steam outlet temperatures. 

5. Provide a framework for future studies involving a 
development of a three dimensional model of a WHRU and a 
detailed, complete COGAS system cycle analysis and optimization. 
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II. MODEL DESCRIPTION 



A. OVERVIEW 

The waste heat recovery unit chosen was a once-through 
counter-cross-flow heat exchanger. A computer program was 
written to model this type heat exchanger for analysis of 
the effect of uniform and non-uniform inlet gas mass distri- 
butions. Variation of other parameters, other than the gas 
distribution, provides a thorough understanding of the oper- 
ation of the heat exchanger in a multitude of conditions. 

The model was written for incorporation into a RACER system 
analysis, but it is not limited to that application. 

Figure 1 shows a three dimensional representation of a 
typical layout of a marine gas turbine waste heat recovery 
unit. The unit is header ed on the fluid inlet and exit and 
designed to be installed in the gas turbine exhaust gas 
duct. The exhaust gas entering the bottom of the WHRU pro- 
vides the energy for transfer to the fluid-side of the system. 
The fluid (feedwater) enters at the upper header of the WHRU 
and flows alternately back and forth through the unit until 
exiting as steam at the lower header. The U-bends shown in 
Figure 1 at the ends of the tube passes allow each tube to 
enter and exit the heat exchanger as a continuous path for 
fluid flow. As a tube makes one horizontal traverse across 
the WHRU it is referred to as a pass. Therefore, the number 
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of times a tube passes through the gas flow defines the 
vertical size of the heat exchanger as a 15, 20, 25, etc., 
pass heat exchanger. 

A gas flow distribution is depicted as entering the WHRU 
in Figure 1. The gas flow is shown oriented along the txibe 
pass as a non-uniform distribution. Additionally, the gas 
flow may be non-uniform in the direction into the heat 
exchanger, i.e., along the axis of the headers as shown in 
Figure 1. This two dimensional gas flow distribution may be 
non-uniform in both mass and temperature. The nature of the 
heat exchanger, by virtue of the cross flow of the fluid, 
decrees that, even if the mass and temperature distributions 
entering the heat exchanger were uniform in both directions, 
these distributions can become distorted and non-uniform as 
the gas passes through the heat exchanger in the third 
dimension. This distortion can be caused by the variation 
of the heat transfer characteristics along the tube. The 
problem of analysis becomes more complicated if cross mixing 
of the gas is considered. 

The model of Combs [Ref. 3J treated the gas distribution 
as a uniform and constant value in one direction. Combs 
assumed the mass and temperature of the gas along the axis 
of the tube was everywhere the same for a complete pass. 

The model yielded results which were close to the results 
attained by other studies and gave a reasonable estimate of 
the overall heat transfer characteristics of the heat 
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exchanger. The model, however, could not predict the effect 
of non-uniformity of the gas distribution on the performance 
of the WHRU. This study discretizes the dimension along the 
axis of the tube pass and enables the designer to evaluate 
the effects of variations in mass flow rate and temperature. 
It is assumed, however, that the distributions along the 
axis of the headers are all the same as the outer distribu- 
tion. The concept of discretization and "nodal" arrangem.ent 
is shown in Figure 2. The representation shows a matrix 
arrangement where the number of rows corresponds to the 
number of passes in the heat exchanger, and the number of 
columns corresponds to the number of gas paths in the 
discretization. The number of passes is a user defined input 
to the model but the number of gas paths is fixed in the 
model at ten. 

Use of the model requires that the user provide the 
following conditions. First, the number of passes that a 
particular water/steam tube is to make through the heat 
exchanger must be specified. This parameter will physically 
size the heat exchanger in the vertical direction. The mass 
flow rate of the gas turbine exhaust gas, (m ), for each of 
ten exhaust gas paths through the heat exchanger must be 
entered. For a uniform distribution, this is simply the 
total exhaust gas mass flow divided by ten. A simplifying 
assumption of the model is that there is no cross mining of 
flow as the gas passes through the heat exchanger; therefore, 
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the mass entering a particular gas path is constant through 
the heat exchanger. 

Next, the temperature of the gas entering each gas path 
is entered. For this study the temperature distribution of 
the gas was considered uniform for all runs; however, the 
model is not constrained to uniform temperature distributions. 

The gas was assumed to be at atmospheric pressure 
entering the heat exchanger and gas side pressure drops were 
not considered in the calculations. 

The user next enters the feedwater inlet temperature, 
the steam outlet temperature and the fluid-side operating 
pressure. These parameters are used in the initial calcula- 
tions to determine a first "guess" mass flow rate on the 
fluid-side of the heat exchanger. The minimum allowed 
exhaust gas temperature is entered. This value is a user 
desired constraint and is specified as the minimum acceptable 
temperature the user can accept to prohibit hot side corro- 
sion. The model is not constrained to maintain the exhaust 
gas temperature above the user specified temperature. 

The number of feet in length, per node, is entered. 

This parameter is based on the physical size of the heat 
exchanger being analyzed in the horizontal direction. For 
all runs in this study, an assumed tube row length of ten 
feet divided by ten gas paths yielded a tube length per node 
of one foot. There is no constraint on the length parameter, 
within reasonable limits. 
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The heat exchanger model was written as a two dimensional 
analysis, i.e., all parameters in the third dimension were 
assumed uniform. The next parameter entered is the number 
of tubes per row, which is used to size the heat exchanger 
in the third dimension. Figure 1 shows a three dimensional 
representation of the heat exchanger design. Each particular 
tube enters at the top of the heat exchanger and passes 
alternately back and forth below itself until exiting the 
heat exchanger at the bottom. The tubes in each subsequent 
plane, behind the first, enter and exit in the same manner, 
but alternate slightly up and down to give a staggered tube 
arrangement. The nodes are numbered in standard matrix format 
with the upper left being node (1,1) , and counting across 
row one as (1,2), (1,3), etc., to (1,10); then node (2,1) 

is directly below node (1,1) , and so on, to the bottom of 
the matrix. Figure 2 shows the nodal arrangement of the 
heat exchanger. The fluid may enter from either side of 
the heat exchanger. If the heat exchanger has an even number 
of passes, the fluid will enter on the same side as the 
steam exit, and on the opposite side for an odd number of 
passes. Each tube is continuous from entry to exit. 

The final user entry is the geometric scaling factor 
which is used in the geometry portion of the model to scale 
all the physical dimensions of the tube, with the exception 
of the length, up or down. 
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Two additional simplifications were assumed in the model. 
First, the fluid side pressure drop was not considered. 
Lastly, even though the tube passes were discretized to 
determine the temperature distributions, the tube length 
within a node still has a significant length for a tempera- 
ture gradient. For analysis, the gas temperature at any 
node was assumed to be uniform across the length of the node 
at both the entry and exit. 



B. INITIAL CALCULATIONS 

The WHRU is essentially analyzed in three sections; 
superheating, boiling and heating. Given an inlet gas distri- 
bution, water/steam inlet and outlet temperatures, and an 
operating pressure, there is a unique fluid-side mass flow 
rate which will satisfy all boundary conditions. The model 
attempts, through iteration, to determine the resultant 
fluid-side mass flow rate. The analysis is also greatly 
affected by the number of passes (i.e., the physical size) 
of the heat exchanger. There are clearly combinations of 
parameters for which no solution is possible, or the solu- 
tion is trivial in that the mass flow rate of the fluid is 
so small as to be meaningless. 

For calculations, variable assignment is as follows; 

(1) m (m,n) = gas flow rate through gas path n at tube 

^ pass m; 

(2) T (m,n) = gas temperature into node (m,n) . The 

^ gas temperature out of node (ra,n) will 

be the gas temperature into node (m-l,n); 



24 



(3) T^(m,n) = fluid temperature into node (m,n) . The 

fluid temperature out node (m,n) will 
equal the fluid temperature into node 
(m,n+l) if the fluid is flowing to the 
right (see Fig. 3) , or the fluid 
temperature out of node (m,n) will equal 
the fluid inlet temperature of node 
(m,n-l) if the fluid is flowing to the 
left . 

At the end of a fluid pass the fluid 
out of a node will equal the fluid 
temperature into the node directly below, 
i.e., T^(m,n) = T^(m+l,n); 

(4) m^ = fluid-side mass flow rate 

Figures 3 and 4 are representations of typical nodes 
depicting the above explanation of fluid and gas flow. 

After all entries are made the program initializes by 
summing the individual mg(NR,n) elements to find total' 
then the T (NR,n) elements are averaged. The T (NR,n) ele- 
ments are averaged by multiplying each element by the corres- 
ponding fraction of the total gas mass that is passing through 
each gas path. These "weighted" temperatures are summed and 
the result is divided by the number of gas paths. The result 
is then the average temperature of the gas into the heat 
exchanger. An average gas temperature is determined using 
the user specified minimum gas temperature out and the aver- 
age gas temperature into the heat exchanger. Then: 



Q = in^4-iC(T .-T .) 

g total pg g avg in g avg out 



where C is evaluated at the heat exchanger gas bulk 
pg 

temperature. Now, a first "guess" fluid mass flow rate is 
calculated from an energy balance 
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Q 



m. 



f out f 



in 



where . and , are the fluid inlet and outlet 

r in r out 

enthalpies at the user specified feedwater inlet temperature 
and steam outlet temperature. 

C . GEOMETRY 

The geometry portion of the model is nearly identical to 
that used by Combs [Ref. 3] . A fin-tube with helically 
curved extended surfaces on circular tubes is assumed for 
this model. The fins are segmented and the tubes are 
configured in banks of one row each with the rows staggered. 
The base tube surface is taken from Weirman, Taborek and 
Marner [Ref. 5] . The description of the finned tubes is as 
follows : 



d. 

1 


= 


tube inside diameter 




zz 


1.86 


in . 


d 

o 


= 


tube outside diameter 




= 


2.00 


in . 


d 

r 


= 


fin root diameter 




= 


2.00 


in. 


^f 


= 


fins per inch 




= 


5.94 




£ 


= 


fin height 




= 


1.015 


in . 


£ 

c 


= 


length of cut from fin 


tip 


= 


0.82 


in. 




= 


fin outside diameter 




= 


4.03 


in . 






fin thickness 




= 


0.048 


in. 


w 

s 


= 


fin segment width 




= 


0.17 


in. 


N 

.Q 


= 


number of segments in 


360 degrees 


= 


38 
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The finned tube configuration is shown in Figure 5. The 
tube length, L, and the number of tubes per row, are 

the user specified entries. The tube layout is shown in 
Figure 6. The center to center tube spacing in the trans- 
verse direction is 4.50 inches. The equilateral, staggered 
tube arrangement used in this model leads to a spacing 
normal to the gas flow, s^, of 4.50 inches and a spacing 
parallel to the gas flow, s^, of 3.90 inches, as shown in 
Figure 6. The heat exchanger height is established by the 
number of WHRU passes (rows) and the tube spacing parallel 
to the gas flow, s . 

In order to establish the minimum gas flow cross-sectional 
area the total "blocked" frontal area, of the heat 

exchanger must be calculated from 



= I- do ^ «t/r " ^ 



and the minimum gas flow area is 



A . = A. - A. 

min f D 

where is the total frontal area of the heat exchanger 
exposed to gas flow. The total inside area available for 
heat transfer per pass is 




7T d . 
1 



L N 



t/r • 
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The gas-side surface area available per pass for heat 
transfer is the sum of the fin surface area and the bare 
tube surface area per pass. The fin surface area per tube 
is calculated from 



A.. 

fin 



N^L[N {2 I w + 2 a + w t ) + ^(d^. - d^)3 

fs cs fc sf 2 fb o 



where d_, = d- - 2£ . The bare tube area per tube is 

fb f c ^ 



A, , = rd L-Trd t^N^L. 

bt o off 

Therefora, the total outside area per pass available for 
heat transfer is 



A 



op 



= Vp**fin ^ 



The cross sectional area for fluid flow is calculated from 



A 



ff 




N 



t/p 



With the mass flow rates, terminal temperatures and heat 
exchanger geometry established, the remainder of the model 
may be solved by nodal analysis for the unique fluid mass 
flow rate and the location of the beginning of the boiling 
and superheating regimes. 
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D. SUPERHEATER NODAL ANALYSIS 



In the heat exchanger there is only one node where two 
conditions are initially known, and that is at the node 
where the steam is exiting. At this node, both the gas inlet 
mass flow rate and temperature and the fluid outlet tempera- 
ture have been specified by the user. Using the initial 
values of the fluid mass flow rate from the estimated overall 
energy balance and the desired superheater steam outlet 
temperature, the heat exchanger is analyzed in the reverse 
direction through the superheater, through the boiler and 
through the heater to arrive at the feedwater inlet. For 
the given geometry and flow conditions, the model will calcu- 
late the water inlet temperature consistent with the assumed 
estimated mass flow rate. If the temperature arrived at by 
analysis as the inlet to the first node does not match the 
user specified feedwater inlet temperature, then the fluid- 
side mass flow rate is adjusted up or down and the entire 
process is begun again. 

To begin, at the final node the gas inlet temperature 
and mass flow rate are known. Also, the fluid outlet tempera- 
ture and a first "guess" fluid mass flow rate are known. 

It is then assumed that the bulk (average) fluid and gas 
temperatures for the node are equal to these known tempera- 
tures. All fluid and gas properties for a node are evaluated 
at the bulk temperatures. 

The gas-side Reynolds number for the node is calculated 
using the bulk temperature. With the Reynolds number, a 
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j -factor is obtained from a polynomial fit to the data for 
ttibe layout number 5 in Ref. 5 (Fig. 6) . By definition, 
the j-factor is related to the heat transfer coefficient, 

hg- by 



j 




By introducing 



St = 



Nu 

Re Pr ' 



the previous expression can be written as 



: = 



Nu „ 2/3 
Pr ^ 



Re Pr 

g g 



Nu 



Re ' 

g g 



and 



1/3 

Nu = 1 Re Pr 

J g g 



where 



Nu 



h 

g 

K 



d 

o 



g 



and K = the thermal conductivity of the gas. Therefore, a 

g 

relationship may be written for the heat transfer coefficient 
as follows; 
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= D 



K 

d 

o 



Re Pr 



1/3 



In the superheat region, the water-side heat transfer 
coefficient is calculated by using the Dittus-Boelter 
correlation 



h 



f 



0.023 




Pr 



0.4 

f 



where all properties are obtained at the bulk temperature of 
the steam in the node. Using the tube wall resistance 
together with h^ and h^, the overall heat transfer coeffi- 
cient for the node can be written as 



U . (m,n) 
oi 



T. £n(d /d. ) 
ip o 1 

2tt K , L 

w t/p 



- + 



A. 

J-P : 

A n. 
op t 



h 



where ; 



n, = 1 - (1-n^) 



A.. 

fin 



f' A 



op 



The fin efficiency, n^, is calculated from the expression 



tanh ML 
ML 



where 



M 



/h P/KA 

g 
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and, if the fin is approximated by a set of rectangular 
strips extending from the tube wall, 



L = 




Now, the cross-sectional area. A, of a rectangular strip 
may be written as 



A 



w 



s 



t 



f 



and the perimeter is 



P = 2 (w + t-) . 

s f 

So, using known geometric parameters and the thermal conduc- 
tivity of the fin metal, the parameter ML may be expressed 
as 



ML = C /h 

g 



where 



C 







Now, 



= [1 - (1 - 



tanh C/h~ 



C/h 



fin, 

A ^ 
op 
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The pass effectiveness is calculated using the effective- 
ness-NTU method. The correlation for the effectiveness for 
a single pass, cross flow heat exchanger with both fluids 
unmixed is taken from Incropera and DeWitt [Ref. 6] 



£ =1- exp[(pj-) (NTU)®’^^{exp[-C (NTU) °’"^^]-l}] 

p C r 



where : 



min 



max 



NTU = 



U . (m,n) A. 
Oi ' ip 

C . 
min 



C . is the lesser of and C where 
min f g 



C, = C ^ m^ 
f pf f 



and 



C 

g 



C m (m,n) 

pg g 



and C is the greater of the two. All fluid and gas 
max ^ ^ 

properties are again evaluated at the bulk temperature. 

Two equations must be satisfied to find the solution to 
a node. These are the rate equation 



Q 



£ C . (T . 
min g in 



■f in 



) 
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and the energy equation 



Q = C ^ (T , , - . 

f pf f out f in 



In the two equations the only unknown is Equating 

the two and solving for yields 



■f in 



^ ^min '^g in ~ ^pf '^f put 

£ C . - ift- C ^ 

min f pf 



Now, Q for the node is 



Q = m,: C ^ (T^ , - . ) 

f pf f out f in' 



and 



T 

g out 



T 

g in 




Next, new bulk temperatures for the node are calculated as 



T , , , = (T . + T /2 

g bulk g in g out 



and 



'^f bulk 



(T 



f out 



T. • ) /2 

f in ' 



An average outside tube wall temperature is calculated from 
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T 

wo 



'^g bulk 






■fB 



) 



where total inside/outside areas for the node 

and are the new node bulk temperatures . The expres- 
sion for is derived from the following formulation for 

heat transfer in the node 



fB _ gB wo 
Z, th o 



which reduces to 



T 

wo 




R 

— 2_(t 

I \h V 



■fB 



) 



where 



R 



n . h A. 

't g to 



and 



I R-i 



U . A. . 
01 ti 



The gas-side film temperature is 



T r, + T 
= qB wo 

gf 2 
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This gas-side film temperature is now introduced at the 
beginning of the calculations for the node as the new gas 
bulk temperature. 

The assumed bulk temperatures, at the beginning of the 
calculations, were obviously too high since the gas inlet 
and fluid outlet temperatures were used. This assxamption, 
however, allowed a solution for the node to be found. The 
model stores the result of this first pass and returns to 
the beginning utilizing the new bulk temperatures to calcu- 
late a new solution to the node. At the end of the second 
pass through, the result of the calculation for T^ is 

compared to the result obtained the first time. If T^ 
has changed, then the new result is stored and the process 
is repeated. When no longer changes, the process has 

converged to a good solution for the node. 

A check is made to determine if the fluid inlet tempera- 
ture is less than or equal to the saturation temperature for 
the pressure specified. If the temperature is higher than 
the saturation temperature, the model proceeds to the next 
node of the superheater and repeats the entire process. 

If the fluid temperature equals the saturation temperature, 
superheating has theoretically begun at the interface between 
this and the final node of the boiling portion of the model. 

If this condition exists, the model proceeds to the boiling 
nodal analysis. 
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As mentioned previously, the model works through the 
heat exchanger in the direction opposite to fluid flow. 



When it is determined that the fluid temperature into a node 
is below the saturation temperature, the node is "split" to 
determine the physical location in the node where super- 
heating begins . 

To split the node the m^ quantity into the node and the 
physical length, L, are halved. The geometry is recomputed 
as if the node were only half as large. The solution 
process is reinitiated and, when a solution is achieved, the 
saturation temperature check is redone. If the temperature 
into the node is still below the saturation temperature, m^ 
and the length are halved again and the process is repeated. 
This procedure is continued until the temperature into the 
node is found to equal the saturation temperature or be 
above it. If the temperature is above the saturation temp- 
erature, then the m^ and L parameters are changed to a value 
halfway between the current and previous values, and the 
process is continued until the temperature into the partial 
node exactly equals the saturation temperature. The values 

for Q, T f T_ . and U . are stored to be used later to 

^ oxi'c X in oi 

develop weighted averages with the boiling portion parameters, 

The remainder of the node is solved in the boiling 
portion of the model . 



E. BOILING NODAL ANALYSIS 

Nodes in the boiling section are solved by the same 
iterative scheme developed in the superheating section. The 
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correlations utilized to determine the solution are differ- 
ent since flow in the boiling regime will involve two phases 
on the fluid-side. In two phase flow the inside heat trans- 
fer coefficient will change with quality, x, and the heat 
flux, q”. The correlation selected for this model is the 
same as that used by Combs [Ref. 3] which was recommended 
by Tong [Ref. 7] for both nucleate boiling and forced convec 
tion. The analysis presented by Combs is presented here. 

The equation for the two-phase flow heat transfer coeffi 
cient was given by Schrock and Grossman in Tong [Ref. 7] as 



h 



TPf 



= B[ 



G h 



fg 






3 

The constants are given by Wright [Ref. 7] ; B = 6.70 x 10 , 

-4 

A = 3.5 X 10 , and n = 0.66. The heat transfer coefficient, 

assuming a totally liquid flow is 



K d. G(l-x) 

= 0.023 ^l— 

1 



0-8 



K, 



-] 



and the Martinelli parameter is 



tt 



[ ^ ] 
^1-x^ 



0.9 P, 0.5 y, 0.1 

[— ] [— ] 

P Mn 

V a 



The limits of the data for this correlation are: 
quality: 0.05 — 0.57 

heat flux: 6 . 0 x 10*^ — 1 . 4 5 x 10 ^ BTU/hr-ft^, 
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In this model, heat flux in the boiling section will be at 

4 2 

levels considerably below 6 x 10 BTU/hr-ft and the full 
range of quality must be considered. Therefore, in order to 
predict the performance of the Schrock and Grossman correla- 
tion at lower heat flux and quality above 0.57, a plot was 
made (Fig. 7) , of the ratio ^TPf'^^L quality where h^ is 

the heat transfer coefficient from the Dittus-Boelter corre- 
lation for water at saturated liquid conditions. The corre- 
lation of Dengler and Addoms , as given in Tong [Ref. 7], 
indicates that low heat flux and high water-vapor mixture 
velocity favor the forced convection mechanism. After a 
review of the information available in Tong, it is not 
entirely clear whether the forced convection or the nucleate 
boiling mechanism is dominant for the relatively low flow 
velocities and low heat flux of this model. In any case, 
the Schrock and Grossman correlation applied at low heat 
flux represents a type of compromise between the entirely 
forced convection assumptions of Dengler and Addoms and the 
"mixed” assumption of Schrock and Grossman. 

A simplification assumed in the model is that the 
quality in the node is the average quality for the entire 
length of the node. Also, to expedite the numerical proce- 
dure, it is assvimed that the node effectiveness is 0.25, 
initially. The effectiveness rarely attains a value higher 
than this for this type heat exchanger, so the assumption is 
valid. The value of the effectiveness is adjusted to the 
"true" value in the iterative process. In a process involving 
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a change of phase the fluid is isothermal. In actuality, 
there is a slight temperature change through the boiling 
regime but for purposes of analysis this change is neglected 
For this reason the fluid bulk temperature in any node in 
the boiling portion is equal to the fluid saturation tempera 
ture . Therefore, to determine the amount of energy absorbed 
by the fluid in a node, the change in fluid enthalpy between 
the inlet and outlet is calculated. 

In the first node in the boiler, which is the node imme- 
diately before the superheater or a partial node if the node 
was split, the gas bulk temperature is again set equal to 
the gas inlet temperature. The gas-side heat transfer coeff 
cient is calculated using the same procedure as delineated 
in the superheater nodal analysis. Next, a Q for the node 
is calculated from 



Q 



e 



m 



g 



C (T 
pg g in 



■f in 



) 



where the effectiveness, s, is initially assumed to be 0.25 
and T^ is equal to the fluid saturation temperature. A 
fluid inlet enthalpy is derived from 



h 



f in 



satv 






where h^ satv ~ enthalpy of the fluid at saturated vapor. 
The fluid bulk enthalpy is 
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fB 



(h 



f in 



h . ^ )/2 

f satv ' 



The average quality for the node is calculated from 



X 



avg 



^fB ~ ^f satl 
^f satv ^f satl 



where h 

liquid 

from 



f satl 
condition. 



the enthalpy of the fluid at a saturated 
The Martinelli parameter is determined 



X 



= [ 



tt 



X 0.9p„0.5y 0.1 

avg , 

'■m/ 

avg V X 



where p^/p^ are the liquid/vapor densities at the saturation 
temperature and y^/y^ are the liquid/vapor viscosities at 
the saturation temperature. Now the heat transfer coefficient 
for totally liquid flow is calculated 



h, 



= 0.023 



d. 

1 



d.G(l-x ) 0.8 C 0 y„ 0, 
[ 1 avg j 






K, 



From this the heat transfer coefficient is calculated for 
two-phase flow 



h, 



TPf 



fg tt 



where h^. = h- . - h- , . and the constants A, B and n 

fg f satv f satl 

have the aforementioned values and G is the mass flow rate 
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per unit area. Finally, the node overall heat transfer 
coefficient is calculated from 



U . 
oi 



h 



TPf 



A. 5,n (d /d . ) 
ip o 1 

2tt K , L 

w t/p 



A. , 

A h 
op g 't 



Whenever two phase flow is involved, the heat capacity of 
the fluid is assumed infinite. For this reason, within 
the boiling regime, will always equal the heat capacity 

of the hot side, so. 



C . 
min 



m 



g 



C 

pg 



and the Number of Transfer Units, NTU, is 



NTU 



U . A. 
oi ip 

C . 
min 



The node effectiveness is 



£ = 1 - exp ( -NTU) 



Now, a revised Q for the node is calculated from 



Q = £ C . (T . - T. . ) . 

min g in f in 



From this Q a revised gas outlet temperature is calculated 
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T 

g out 



T 

g in 



‘0/=min> 



With this gas temperature out, a new film temperature is 
calculated in the same manner as for the superheating section. 

The gas outlet temperature is stored and the film temp- 
erature is used to repeat the process for the node. At the 
end of the second time through, the gas outlet temperature 
derived is compared with the result from the first. If 
equality exists, a solution for the node has been found and 
the model proceeds to the next node. If inequality exists, 
the result of the current iteration is stored for comparison 
with the results of the following iteration. 

When a final solution is attained, the enthalpy of the 
fluid at the inlet to the node is checked against the 
enthalpy of saturated water at the heat exchanger pressure. 

If the inlet enthalpy is greater than or equal to the enthalpy 
of saturated water, the model proceeds to the next node and 
begins the solution process again. If the inlet enthalpy 
is less than the enthalpy of saturated water then the model 
proceeds to a node splitting routine similar to that used 
when transitting from the superheating regime to the boiling 
regime. 

The node splitting is again accomplished by halving the 
mass flow rate on the gas-side and halving the physical 
length, L, of the node. The remaining geometric parameters 
are recomputed and the solution process is reinitiated. The 



43 



mass and length adjustment procedure is continued until the 

inlet enthalpy is equal to the enthalpy of saturated water. 

When this is achieved, the values of T , U . and Q are 

g out oi 

stored to be recombined as weighted averages with the heating 
portion of the node to give the overall nodal values. 

The remainder of the node is solved in the heating portion 
of the model . 

F. HEATING NODAL ANALYSIS 

The procedure for analyzing a node in the heating portion 
of the model is identical to that used in the superheating 
portion. Each node is analyzed in an iterative process, and 
it makes no difference if the first node was a split node. 

The first node in the heating section will have a gas inlet 
temperature equal to the gas outlet temperature of the node 
directly below. The fluid outlet temperature will equal the 
fluid saturation temperature for the specified pressure. 

An outline of the calculations is as follows; 

1. Calculate the gas-side heat transfer coefficient 
using the node gas inlet temperature as the gas bulk 
temperature . 

2. Calculate the fluid-side heat transfer coefficient 
initially using the fluid saturation temperature as 
the fluid bulk temperature. 

3. Calculate the overall heat transfer coefficient using 
the ef f ectiveness-NTU method. 
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4, 



Determine the node fluid inlet temperature from 



•f in 



eC. T. -m_C_T^ , 

min q in f pf f out 

e C . - m^ C 

min f pf 



5. Calculate a first "guess” for Q for the node from 



Q = m^ C .(T. . - T. . ) 

f pf f out f in' 



6. Calculate the node gas outlet temperature from 



T . = T . - (Q/m C ) . 

g out g in g pg 



7. Calculate T , ,, = (T . + T /2 

g bulk g in g out ' 



8. Calculate a node average film temperature from 



■gf 



= (T 



gB 



T )/2 
wo 



9. Store the result of the first pass through the model 
and iterate until procedure converges to solution, 
utilizing gas outlet temperature for comparison at 
each pass. 

All properties of the gas and fluid are evaluated at the 
respective bulk temperatures for the node. The process is 
repeated node by node until the node at the fluid inlet to 
the heat exchanger is completed. The final node is either 
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node (1,1) for an odd number of passes, or node (1,10) for 
an even number of passes in the heat exchanger. 

G. HEAT EXCHANGER SOLUTION PROCEDURE 

When the procedure of the heating section has arrived at 
a solution for the final node of the heat exchanger, the 
fluid inlet temperature for the node is compared to the user 
specified feedwater inlet temperature. If equality exists, 
the model proceeds to an output routine where the matrices 
of water/steam temperatures, exhaust gas temperatures, 
overall heat transfer coefficients and individual heat 
transfer per node are printed. Additionally, the fluid-side 
mass flow rate, minimum pinch point temperature difference, 
total Q for the heat exchanger and isentropic steam turbine 
power are calculated and printed. 

The pinch point temperature difference is defined in 
this model as the minimum difference in the heat exchanger 
between the gas inlet temperature and the fluid inlet temp- 
erature to any node. Consideration of the effect of a pinch 
point below approximately 25 deg. F is cause for concern 
only if the low pinch point is coupled with a low gas temp- 
erature. This combination is non-prohibitive to hot side 
corrosion. If the final node fluid inlet temperature is not 
equal to the user specified feedwater inlet temperature, then 
the fluid mass flow rate that was "guessed" at the beginning 
is adjusted and the entire model is run again. If the solution 
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temperature is too high the "guessed" fluid mass flow rate 
is lowered. With less mass going through the heat exchanger 
the flow velocity is decreased. A lower flow velocity through 
the heat exchanger enables more energy from the gas-side to 
be transferred to the fluid-side. A larger Q will increase 
the temperature differential between the fluid inlet and 
outlet which will yield a lower solution fluid inlet temp- 
erature the second time through the model. The process is 
opposite, i.e., mass is added, if the fluid inlet temperature 
to the first node is found to be too low. 

The model will continue to adjust the mass after each 
iteration in smaller amounts until an exact match is made to 
the user specified feedwater inlet temperature. When equality 
is achieved, the model proceeds to the output routine where 
the above mentioned parameters are printed. 

Included, as Appendix A, are sample outputs of the model 
for the uniform distribution, distribution 1 and distribu- 
tion 5 for all four of the selected gas mass flow rates. 

The run numbers at the top of the first page of each output 
were utilized for accounting purposes. The uniform distri- 
bution was used for runs 1, 10, 19 and 28; distribution 1 was 
used for runs 2, 11, 20 and 29; and distribution 5 was used 
for runs 6, 15, 24 and 33. The first page of the output shows 
the inter-active questions asked by the model and the user 
responses. The following pages show the matrices developed 
by the model in the solution process. Comparison of the 




mmmmrn ^ ^ 



i 
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matrices for the uniform and non-uniform distributions 
demonstrates the serious impact of the non- uniformity on the 
performance. The non-uniform distribution profiles can be 
seen as projected profiles on the results of each of the 
matrices. The final page of the output is a listing of 
additional parameters which are calculated by the model includ- 
ing the fluid mass flow rate, Q based on the individual 
matrices, location of the beginning of boiling and super- 
heating within the WHRU and the isentropic horsepower. 
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III. DISCUSSION OF RESULTS 



A. INITIAL TESTING 

The foregoing model was coded in FORTRAN IV and applied 
in a computer simulation on an IBM- 3 033 computer. 

To properly test the model and develop an understanding 
of the operation of the WHRU several runs were made utilizing 
a uniform gas distribution with a total mass flow rate and 
temperature corresponding to a gas turbine horsepower rating 
for approximately a 20 kt. speed. The conditions presupposed 
usage of a General Electric LM2500 Gas Turbine Engine, and 
gas mass and temperature data was taken from [Ref. 8] . Initially, 
the model was run with attention to pinch point temperatures. 

The number of passes was varied to find the largest heat ex- 
changer which could be used and still maintain a pinch point 
in excess of 25 degrees fahrenheit. Under conditions of non- 
uniform flow distributions the model definition of pinch point 
is misleading since small gas mass flow rates through a 
particular gas path may yield a nearly zero pinch point when 
the overall performance of the heat exchanger is not affected. 
Inspection of the model output for the non-uniform distribu- 
tions will demonstrate this result. The preliminary runs 
allowed selection of a 23 pass heat exchanger with a 29.4 deg. 

F pinch point. Therefore, the 23 pass heat exchanger became 
the "design.” The performance data tabulated for the LM2500 
engine [Ref. 8] was utilized to select four gas mass flow 
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rates and inlet temperatures which correspond closely to 
data points on a cubic power curve. The four combinations 
of mass and temperature chosen were; 



m = 

g 


121. 


8 Ibm/s, 


T . 
g in 


= 844 


° F, 


SHP = 


^ 17500, 


NPT = 3000 rpm 


m = 

g 


99.1 


Ibm/s , 


T . 
g in 


= 755° 


F, 


SHP = 


10000, 


NPT = 2400 rpm 


m = 

g 


77.4 


Ibm/s , 


T . 
g in 


= 702° 


F, 


SHP = 


5000, 


NPT = 1800 rpm 


m = 

g 


66.7 


Ibm/s , 


T . 
g in 


= 694° 


F, 


SHP = 


3000, 


NPT = 1200 rpm 



Each of these was run in the model as uniform and various 
non-uniform distributions. 

To develop the non-uniform distributions, a normalized 
flow distribution table was constructed and is presented as 
Table I. This table shows the fraction of the total mass 
flow rate which is passing through each individual gas path 
for each of nine different flow distributions. The flow 
distributions were arbitrarily derived in an attempt to 
develop a spectrum of distributions that may enter the heat 
exchanger. A graphical display of each of the flow distribu- 
tions is shown in Figures 8A-8I. The total mass flow rate 
used was multiplied by each entry in the normalized flow 
distribution table to arrive at the mass flow rate through 
each of the ten gas paths for any of the nine distributions. 
This procedure allowed a comparison of the effects of non- 
uniform distributions on the performance of the heat ex- 
changer. All nine distributions were run for each of the 
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four chosen mass and temperature combinations. The actual 
mass flow distributions used are shown in Tables II-A--II-D. 
The combination of nine distributions and four different 
total gas flow rates yielded 36 runs. For the runs, the 
remaining parameters were: 



a) 


feedwater inlet 


temperature 


= 


200 


deg . 


F 


b) 


desired steam temperature 


= 


650 


deg . 


F 


c) 


steam operating 


pressure 


= 


400 


psia 




d) 


minimum allowed 


gas outlet temperature 


= 


300 


deg. 


F 


e) 


node length 




= 


1.0 


ft. 




f) 


number of tubes 


per row 


= 


32 






g) 


geometric scaling factor 


= 


1.0 







The model calculates the isentropic horsepower that would 
be developed by an ideal steam turbine. This power can also 
be viewed as the maximum power available from the developed 
superheated steam exhausting to a fixed condenser pressure. 

The method of calculation is to determine the enthalpy of the 
steam at the outlet pressure and temperature and then sub- 
tract the enthalpy of the steam after isentropic expansion 
to a condenser pressure of 2 psia. The difference in enthalpy 
was multiplied by the fluid mass flow rate, with the appro- 
priate conversions, to arrive at a horsepower value. The 
horsepower calculated for each run is tabulated in Table III. 
The table lists the results of each run of the model for all 
36 runs for each of 5 different steam outlet temperatures. 
Table IV is a similar table showing the solution values of 
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the fluid mass flow rates derived in each run. The two 
tables demonstrate the variation of fluid mass flow rate and 
power which is attained by varying the desired steam outlet 
temperature; and the effect of the non-uniform distributions 
on performance. 

B. EFFECT OF NON-UNIFORM DISTRIBUTIONS 

The data from column 1 of Table III was plotted and is 
shown in Figures 9A-9E as power vs. gas flow rate. Figure 
9A shows the comparison of the uniform distribution with 
distributions 1 and 2. Figure 9B shows the comparison of 
the uniform distribution with distributions 3 and 4, etc. 

From the graphs it can be seen that the distribution inflicting 
the most serious impact on the heat exchanger performance is 
distribution 1 (see Figure 8-B) . Distribution 1 is highly 
center concentrated and utilizes the area available for heat 
transfer the least efficiently. For comparison, distribution 
3 is a nearly uniform distribution and it inflicts a relatively 
minor impact on the performance. Distributions 2 and 4 are 
ranked next in order of impact. These are also distributions 
where the mass is centrally concentrated, but to a lesser 
degree than distribution 1. Distributions 5 and 6 are mirror 
images of one another and also have a zero mass flow rate 
through two of the gas paths. A zero mass flow rate indicates 
that the flow through that particular pass has no velocity. 

In a real situation there will be a transient heat transfer 
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from the still mass of gas surrounding the tubes by natural 
convection until the gas temperature reaches the tube wall 
or fluid temperature. In any case, the amount of heat 
transferred in such a node is negligible compared to the 
nodes experiencing a gas flow rate. For the model, a zero 
gas flow rate in a node will be treated as having zero heat 
transfer, i.e., a "dead node," therefore, the Q for the node 
was set to zero, the fluid inlet temperature and fluid outlet 
temperature were set equal and the gas outlet temperature 
was set equal to the gas inlet temperature. In reality, the 
gas temperature in a "dead node" will attain the temperature 
of the fluid in the node but for pinch point calculation pur- 
poses the original gas inlet temperature was maintained. 

Also, since the model assumes no cross-mixing of the gas 
between gas paths as the distribution passes through the 
heat exchanger, this assumption ignores the potential tendency 
of the gas mass distribution to become uniform during passage. 
Distributions 5 and 6 are skewed to the left and right, 
respectively. The combination of the null gas flow paths and 
the skewedness of these distributions caused the second most 
detrimental effect on the heat exchanger. It should be noted, 
however, that whether the distribution was skewed to the right 
or left side had relatively no effect; that is, the results 
for distribution 5 are nearly identical to those for distri- 
bution 6. The same result was true for distributions 7 and 
8, which are also mirror images of each other with left and 
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right skewedness. The effect of distributions 7 and 8 were 
very nearly the same as those for distribution 4. Figure 9E 
is a composite which demonstrates the effect of all the dis- 
tributions. The ranking as to effect yields: 

1. the uniform distribution gives the best performance 

2. distr. 3, nearly uniform 

3. distr. 4, least center concentrated, nearly same as 
distr . ' s 7 and 8 

4. distr. 2 , intermediately center concentrated 

5. distr. 's 5 and 6, skewed with null flow paths 

6. distr. 1, highly center concentrated 

The ranking can be viewed as a measure of the efficiency 
of the distribution in utilizing the area available in the 
heat exchanger for heat transfer. As the amount of distor- 
tion of the flow increases, the quantity of energy trans- 
ferred decreases. 

Figures lOA-lOE are graphical representations of the 
solution fluid mass flow rates vs. the gas mass flow rates 
for each of the distributions; exactly as was done to develop 
Figures 9A-9E. As was to be expected, the results for Figures 
lOA-lOE are identical to those for power vs. gas flow rate, 
since power is directly proportional to the fluid mass flow 
rate. 

After these preliminary results, it was decided to con- 
duct the remaining runs utilizing only the uniform distribu- 
tions and the worst of the non-uniform distributions; i.e.. 
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distribution 1. It was assumed that results for the remain- 
ing distributions would fall between the values of these two. 

C. EFFECT OF HEAT EXCHANGER SIZE 

Next, an analysis was done to determine the effect on 
output power and fluid mass flow rate by varying the number 
of heat exchanger passes. Using the four chosen gas mass 
flow rates, runs were conducted for heat exchangers of 20, 

21, 23 and 25 passes. All runs were for 650 deg. F outlet 
steam at 400 psia operating pressure. Therefore, the runs 
for the 23 pass heat exchanger were the "design" runs and 
the runs for 20, 21 and 25 passes were the "off-design" runs 
provided for comparison. The results of these runs are 
tabulated in Table V and graphically presented in Figures 
llA-llH as power and fluid mass flow rate, respectively, vs. 
the number of heat exchanger passes. Again, because of the 
relationship between power and fluid mass flow rate, the 
graphical results are similar in contour. The results show 
that as the size of the heat exchanger is increased, so will 
the power out. To provide the increased power the mass flow 
rate of the fluid must increase since for the different 
sized heat exchangers the inlet and outlet enthalpies of the 
fluid is not changing. The limit to the increase in size is 
the available energy in the gas. Increasing the heat ex- 
changer size to the point where sufficient energy is trans- 
ferred from the gas to the fluid to attain the minimum 
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acceptable exhaust gas temperature is the maximum size 
possible. 

The curves display an interesting result concerning the 
differences between the uniform distribution and the non- 
uniform distribution. As the size of the heat exchanger 
increases, the two curves are converging with the difference 
in power (or fluid mass flow rate) being less for the larger 
heat exchanger than for the smaller. From the data of Table 
V, it can be seen that, for the 121.8 Ibm/s gas flow rate, 
the power from distribution 1 is 89% of the power attained 
from the uniform distribution for a 20 pass heat exchanger. 
This increases to approximately 91% for the 25 pass heat 
exchanger. The more rapid increase in power for the non- 
uniform distribution is caused by the mass concentration in 
the center gas paths. For a particular size heat exchanger, 
the paths with a small gas mass flow will attain the temper- 
ature of the fluid, usually in the boiling region. From that 
point on, these paths contribute little, or not at all, to 
the total heat transfer. The high concentration of mass 
in the center paths, however, contribute greatly to the total 
heat transfer because of the high velocity of the flow and 
the high heat transfer coefficients. These paths of high 
concentration will contribute until the gas exits the heat 
exchanger. By adding more passes to the heat exchanger 
it is as if the "dead areas" in the low mass flow paths were 
placed in the path of the high mass flows. The effect is a 
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significant increase in the heat transfer. The high heat 
transfer coefficients, mass velocity and temperature differ- 
entials all contribute to a greater percentage increase per 
pass than is achieved by the uniform distribution. As addi- 
tional passes are added, the performance of the non-uniform 
distribution improves, but it will never attain the level of 
the uniform distribution. For two heat exchangers of equal 
size, the uniform distribution gives the most efficient use 
of the available area for heat transfer. 

D. EFFECT OF STEAM OUTLET TEMPERATURE 

Each distribution, uniform and distributions 1-8, was 
run at each of the four gas flow rates for five specified 
steam outlet temperatures. The power and fluid mass flow 
calculated by the model are presented in Tables III and IV. 

The results are also presented graphically in Figures 12-A 
and 12-B for the uniform distribution and distribution 1 
at each of the four gas flow rates. 

Both the tabular and graphical data shows a nearly uni- 
form increase in fluid mass flow rate with the linear de- 
crease in steam outlet temperature. The tabular data for 
power demonstrates a non-linear result, showing initially a 
decrease in power with a decrease in steam outlet temperature, 
and then increasing to attain a power higher than that for 
a 650 deg. F steam outlet temperature. 

The non-linearity of the power is caused by the non- 
linearity of the steam enthalpy as the steam outlet temperature 
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decreases. The combination of the change in enthalpy and 
the change in the fluid mass flow rate results in a mislead- 
ing impression of attaining a higher pov/er at a lower steam 
temperature. The total power does increase but the specific 
power decreases. The increase in fluid mass flow between a 
650 deg. F steam outlet temperature and a 550 deg. F steam 
outlet temperature is approximately 7%. The increase in 
power at the lower steam outlet temperature is less than 0.2%. 
The penalty to utilize a lower steam temperature and main- 
tain an equivalent power output is a water/steam system which 
must be able to circulate 7% more fluid in the same time 
period. This result may have serious impact in the design 
when condenser size, circulating pump size, valves, piping 
and the steam turbine are taken into consideration. The 
steam turbine must also be considered with respect to the 
lower steam temperature having a lower quality through the 
latter stages of the expansion. The effects of overly "wet" 
steam can be disastrous to many turbines. The end result 
is that the sacrifice in power at the higher steam tempera- 
ture is justified by the decrease in fluid flow required 
and higher end point quality. 

E. EFFECT OF HEAT EXCHANGER PRESSURE 

The uniform distribution and distribution 1 were run for 
each of the four gas flow rates with three different speci- 
fied heat exchanger operating pressures. The pressures were 
the "design" pressure of 400 psia and then 500 psia and 600 
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psia. The tabulated results of these runs are presented in 
Table VI, and graphically presented in Figures 13-A and 13-B 
for power and fluid mass flow rate, respectively, vs. heat 
exchanger pressure. Again, because of the direct correla- 
tion between fluid mass flow rate and power, the contours of 
the graphical results are similar. 

The curves show a decrease in power and fluid mass flow 
rate as heat exchanger pressure is increased. As the pres- 
sure is increased by 50% above the design level, there is a 
4% drop in the system fluid mass flow rate and a 7.5% loss 
in power. Clearly, the benefits derived from the decrease 
in system fluid mass flow rate does not justify the sacri- 
fice in system power. With respect to the analysis of the 
previous section, a compromise must be reached to achieve 
the maximum power from the system within the confines of a 
maximum acceptable fluid mass flow rate. 

The graphical results also show the effect of the non- 
uniform distribution for each of the four gas flow rates. 

As the gas flow rate decreases, the penalty for non-uniformity 
has a lesser effect. The percentage loss for the four flow 
rates are approximately 9, 8, 7 and 6 percent, respectively, 
from highest flow rate to lowest. The loss is approaching 
zero as the flow rate approaches zero. The loss in power is 
directly related to the decrease in fluid mass flow rate 
determined in the solution; and the decrease in fluid flow 
rate is a result of the ineffectiveness of the non-uniformity 
on the operation of the heat exchanger. 
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F. EFFECT OF GAS FLOW RATE ON HEAT EXCHANGER PERFORMANCE 

The data from the matrix of water/steam temperatures for 
the four gas flow rates as uniform distributions was plotted 
as a function of tube length. The curves are shown in 
Figures 14-A — 14-D. Figure 14-E is a plot of the water/ 
steam temperatures as a function of tube length for distri- 
bution 1 with a gas mass flow rate of 121.8 Ibm/s . Figure 
14-F is a superposition of Figure 14-E on Figure 14-A. 

Analysis of Figures 14-A through 14-D shows the effect 
of the decrease in gas flow rate on the rate of heat trans- 
fer as the fluid passes through the heat exchanger. The 
higher the gas flow, the steeper is the slope of the curve 
in the superheat region. A steeper slope is indicative of 
a shorter tube length required to attain the desired super- 
heat temperature. The effect, as the gas flow rate is 
decreased, is that more of the heat exchanger is dedicated 
to superheating. The length of the tube required to raise 
the water from a saturated liquid to a saturated vapor 
condition is nearly the same for all gas flow rates with 
only a ''left shift" on the plots caused by the superheat 
conditions. The slope of the curves in the heating region 
is the opposite of that of the superheat. As the gas flow 
rate decreases the slope of the curve in the heating region 
increases . 

Since, for all conditions, the feedwater inlet tempera- 
ture and steam outlet temperature are the same, all ti±>es 
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are continuous and of the same length, the difference in the 
curves is attributable to the rate of heat transfer and the 
total amount of energy transferred. The rate 'equation is 
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The change in gas mass flow rate will obviously decrease the 
total amount of energy transferred since less gas in means 
less energy into the heat exchanger. If the curves were 
plotted in a third dimension with Q, the smallest Q would 
occur with the smallest gas flow rate. The forced equality 
with the rate equation means that, as the gas flow rate 
decreases, the rate of heat transfer must decrease. The 
figures show the total effect on the profile operation of 
the heat exchanger is of little consequence. The similarity 
of the curves is an indication of the "physical operation" 
of the heat exchanger under varying conditions. For the test 
runs, a decrease in gas flow rate was coupled with a decrease 
in gas inlet temperature. The two parameters changing and 
the similarity of the curves emphasizes the physical constancy 
of the design. 
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Figure 14-E is the plot for distribution 1. The "ripple" 
effect in the heating and superheating regions is a result 
of the non- uniformity of the gas mass distribution, but the 
superposition shown in Figure 14-F demons.trates the oblivious- 
ness of the physical properties of the heat exchanger to a 
variation of input. 

G. EFFECT OF NON-UNIFORM MASS DISTRIBUTION ON EXHAUST GAS 

TEMPERATURE 

For all runs of the model the temperature distribution 
of the exhaust gas entering the heat exchanger was assumed 
uniform. To determine the effects of a non-uniform gas mass 
distribution on the exhaust gas temperature profile, a plot 
was made of the gas temperatures across several of the heat 
exchanger passes. The matrix of exhaust gas temperatures 
was used and data was plotted for both a uniform distribu- 
tion and distribution 1. The gas mass flow rate was 121.8 
Ibm/s. Figures 15-A and 15-B show the results. 

Figure 15-A demonstrates that with a uniform mass and 
temperature distribution at the inlet to the heat exchanger, 
the distribution is virtually undistorted by passage through 
the WHRU. This result must be considered ideal with respect 
to the assimption of no cross-mixing of the gas during passage. 
There is a smooth transition of energy from inlet to outlet 
which accounts for the best performance of the heat exchanger 
under this condition. The effect of the non- uniformity is 
clearly seen in Figure 15-B. The temperature distribution is 



62 



rapidly distorted to become almost a replica of the gas flow 
distribution. This is a result of the disparity of heat 
transferred and heat transfer rates from node to node in 
the WHRU. A slight smoothing effect is seen in the boiling 
portion, but distortion is again accentuated in the heating 
section. In a real situation there will be some non- 
uniformity of temperature as it enters the heat exchanger. 

The analysis shows that the condition will be worsened to 
the degree of the mass non- uniformity with the subsequent 
loss of performance in the WHRU. 

The distortion of the exhaust gas temperature distribu- 
tion forces consideration of pinch point. When the tempera- 
ture differential between the exhaust gas and the tube 
becomes less than approximately 25 deg. F coupled with a gas 
temperature less than approximately 300 deg. F, the WHRU may 
experience localized hot-side corrosion due to sulphur products 
in the exhaust gas. It should be noted that this condition 
will exist only near the top of the heat exchanger near the 
exhaust gas exit. If the degree of flow distortion is 
severe and maintained for a long period, the effects can be 
disastrous . 
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IV. CONCLUSIONS 



The waste heat recovery unit is a practical method of 
recouping thermal energy from gas turbine engine exhaust 
gases for fuel economy improvement. In the design and 
installation of a RACER system several factors must be taken 
into consideration. 

1. Mass Flow Distribution. The performance of a waste 
heat recovery unit is adversely affected by a non-uniform gas 
mass distribution; and the degradation is directly propor- 
tional to the degree of non-uniformity. To ensure acceptable 
results, the system should have flow distribution control 
devices installed between the gas turbine and the WHRU. 

2. Operating Pressure. The operating pressure of the 
system affects the power output inversely. A higher pressure 
system requires less fluid to deliver the same power but the 
percent reduction in fluid required does not justify the 
percent loss in power. Since the system is not a main pro- 
pulsion system requiring a large volume of fluid, the model 
suggests using a low pressure system. 

3. Steam Outlet Temperature. Reducing the outlet temp- 
erature required forces an increase in fluid mass flow rate 
to maintain the same power out. High temperature steam, how- 
ever, has a greater wear effect on system components such as 
valves, seals and piping. The model suggests a steam tempera- 
ture as high as practical should be used, within material 
constraints . 
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4. Heat Exchanger Size. Power and fluid mass flow rate 
are directly related to heat exchanger size. The size of 
the WHRU should be as large as the available space and 
weight constraints allow, or as large as necessary to pro- 
duce the power required. Consideration must be given to the 
system as a whole since, as the heat exchanger becomes larger, 
so does the piping, valves, pumps and condenser necessary to 
operate. 

5. Speed. Power developed is a function of the energy 
available in the exhaust gas flow. At higher gas turbine 
speeds more horsepower is available from the WHRU to assist 
the main propulsion system. Utilizing the latest NAVSEA 
speed profile, the heat exchanger should be sized in accordance 
with the ship speed most commonly used. 

All data suggests a WHRU as large as physically possi- 
ble, operated at a low steam pressure, high steam temperature, 
and with flow distribution control devices installed provides 
the best performance. 
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V. RECOMMENDATIONS FOR FURTHER RESEARCH 



A. WHRU MODEL IMPROVEMENT AND EXPANSION 

The control imposed on the WHRU model in this study was 
adjustment of the steam flow rate to maintain a user speci- 
fied steam outlet temperature. This forces the steam system 
to respond to a variation of gas inlet conditions by "throt- 
tling" the feedwater system to maintain the required steam 
temperature . 

This is only one of several types of control systems 
which may be used in a RACER system installation. Two 
possible alternative systems are: one which maintains a 

constant steam flow rate while allowing the steam tempera- 
ture to fluctuate, and one in which pressure is allowed to 
"float" with steam demand. In this second type system, the 
flow resistance of the system is used, and as the fluid mass 
flow rate decreases (steam demand decreases) , the pressure 
required decreases. Knowing the inter-relationship will 
yield a "floating" as well. Investigation into the 

performance of a WHRU with these types of control could lead 
to an improved output. The investigation of all types of 
control should consider the effects on the entire system 
because, the type of control may dictate the type steam 
turbine which must be used. 

To expand the understanding of the effects of non-uniform 
gas distributions on the WHRU performance, the current 
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model could be enlarged to allow a variation of the input in 
an additional dimension. A two dimensional non-uniform input 
will yield a three dimensional picture of the temperature 
distributions as the gas passes through the heat exchanger. 

An analysis of this type would help design the flow distri- 
bution control devices to eliminate severe trouble spots 
contributing to hot side corrosion and to maximize the WHRU 
performance. 

In addition to the expansion to a three dimensional 
analysis, the current model could be improved by including 
the following items. 

1. Fluid-Side Pressure Drop 

The water/steam-side pressure drop may be a signifi- 
cant percentage of the system operating pressure. Calculation 
of the pressure drop and its effect on system performance 
should be included. 

2 . Additional Geometries 

Additional fin-tube configurations should be con- 
sidered for possible improvement of the heat transfer charac- 
teristics and the weight and space requirements of the WHRU. 

3 . Fouling 

Consideration of inside and outside heat transfer 
surface fouling should be given when investigating possible 
fin- tube configurations. 

4 . Alternate Working Fluids 

Working fluids, other than water, should be investi- 
gated for possible enhanced thermodynamic characteristics at 
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a miniituam cost in system maintainability, while maintaining 
safety to operating personnel. 

5 . Non-Uniform Temperature Distributions 

The current model could be utilized to determine the 
effects of non-uniform temperature distributions. The gas 
is at a uniform temperature as it exits the gas turbine, but 
it probably is not uniform at the inlet to the WHRU. Inves- 
tigation into the effects of non-uniformity of both gas mass 
and temperature should be accomplished. 

The model should be incorporated into a complete 
RACER system analysis, including models of the steam turbine, 
condenser and circulating pump. Variations of gas distribu- 
tion input can then be analyzed to determine total system 
effect. Understanding the total system effect and the inter- 
relationship of the variables could then be used to optimize 
the system. 
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NORMALIZED FLOW DISTRIBUTION 
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TABLE II-B 
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1 


2 


3 


A 


5 


6 


7 


a 


Q 


10 


U 


9.91 


9.91 


9. 9l 


9.91 


9.91 


1 9.91 


t 9.91 


i 9.91 


9. 91 


9. Ql 


1 


1 .68 


4 .85 


6. 74 


11.10 


25.17 


1 2S. 1 7 


11.10 


1 6.74 


1 4.85 


1 . 68 


2 


3.06 


4.36 


8. S2 


1 3 . 1 6 


20.4 1 


20.41 


1 13.18 


i 9.52 


1 4.36 


3.08 


J 


7.92 


0.92 


9. 91 


10.90 


1 1 .90 


1 11.90 


10.90 ! 


j 9. 91 ! 


1 8.92 


7. 92 


A 


2.46 


6.44 


1 0. 41 


13.37 


•6.35 j 


i 16.35 


13.87 


1 10.^1 


6.44 


2.48 


5 


9.91 


22.20 


1 9. 03 


16.15 


12.99 


9. Si 1 


1 6.35 1 


1 2-97 i 


1 0.00 


0. 00 


6 


0.00 j 


1 0.00 


2.97 


6.35 


9.5l 


12.99 j 


16. 15 1 


1 19.03 


1 22.20 


9.91 


7 


*2.59 


19.23 


1 S.56 


12.29 


3.52 


7. 23 j 


6.60 1 


6.54 


4.31 


4. 26 


6 


4.26 


6.31 


6.54 


6.50 


7.23 


8*52 


12.29 


15.56 


19.23 


12.59 
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TABLE II-C 



GAS FLOW DISTRIBUTIONS (m 

g 



77.4 LBM/S) 



6*S K*SS FLO* RATE s 77.* LBM/S 

GAS TEMPERAruRE AT MEAT EKCHANCEB INLET = 702 OEG. F 
CASE 5S3 I REF 8 ) SHP « 5000 NPT * I BOO 



OIST* 






gas 


path 














NO* 


t 


2 


3 


4 


5 


6 


7 


a 


9 


10 


U 


7.7 4 


7. 74 


7.74 


7 .74 


7.74 


7.74 


7.74 


7,74 


7.74 


7. 74 


1 


1.31 


3.80 


5.26 


9.67 


19.66 


1 9.66 


8.67 1 


1 5.26 


3.80 


1.31 


2 


2*41 


3. Al 




10.29 


IS. 94 


15.94 


10.29 1 


6.65 


3. 41 


2.41 


3 


6. 1 9 


6.97 


7. 74 


8.51 


9.29 


9. 29 


6.51 


1 7.74 


6.97 


6.19 


4 


1 . 9A 


S.03 


6. 1 3 


10.33 


12.77 

j 


12.77 


10.83 


8. 1 3 


5.03 


1 . 94 


s 


7*74 


17.34 


14.86 


12.61 


10.15 


7.43 


4 .96 


2. 32 


0.00 


0.00 


6 


0 .OO 


0.00 


2. 32 


4.96 


7. *3 1 


10.15 


12.61 


14.86 


1 7.34 


7.74 


7 




1 S.02 


12.15 


9.60 1 


1 6.6S j 


j S.6S 


5. 15 j 


1 5.11 


4.93 


3.33 


3 


3.33 


<• .9J 


1 S.M i 


! 5.15 1 


1 5-65 i 


1 6.65 


9.60 ; 


1 2 . 1 S 


15.02 


9.83 



TABLE II-D 

GAS Flow DISTRIBUTIONS (m =66.7 LB^V'S) 

9 



GAS MASS PL3W RATE = 67.7 LBM/5 

GAS TEMPERATliftE AT m£a T EXCHANGER INLET = 6RA OEG* P 
CASE SA5 (REF 8 ) SHP r 3003 NPT * 1200 



OIST . 






GAS 


P4TM 














NO. 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


ii 


6.77 


1 6.77 j 


1 6.77 


6.77 


6.77 


6.77 


6.77 


6.77 


6.77 


6.77 


1 


1.15 


3.32 j 


1 4.60 


7.58 ; 


1 17.20 


17.20 


7.58 


1 4 .60 


3.32 


1.15 


2 


2.11 


2.96 j 


1 5.82 


9.00 


' 13.9a 1 


1 13.9* 


9.00 ! 


1 5.32 


2.98 


2.11 


3 


S.41 I 


6.10 


1 6. 77 


7,44 


8.13 


6. 1 3 


7.44 


1 6.77 


6. 10 


5.41 


4 


1.70 j 


4.40 


7.11 


9.4 7 


11.17 


11.17 


9. *7 


I 7.11 


A . 40 


1 • 70 


s 


6.77 1 


1 S . 1 7 


1 3. 00 


11.03 


a.aa 


6.50 


4 .3 A 


2.03 


0.00 


0.00 


6 


0.00 1 


1 0.00 1 


i 2.03 


4 .34 


6.50 j 


j 6-98 


11.03 


13.00 1 


1 15.17 


6.77 


7 


d.60 j 


13.14 


10.63 


8.40 


S.82 j 


j 4.94 


4.50 


4.47 


4.31 


2.91 


8 


2.91 


4.31 


4.47 


4.50 


4.9a 1 


S.S2 


8 .AO ; 


1 10.63 


1 3. 1 A j 


6.60 
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TABLE III 



ISENTROPIC POWER FOR VARIOUS GAS FLOW DISTRIBUTIONS, 
GAS FLOW RATES AND STEAM OUTLET TEMPERATURES 





RUN 

NO. 


DISTR. 
NO . 


STEAM TEMPERATURE 

650 625 600 575 550 


GAS 
MASS 
121. 8 
LBM/5 


1 


U 


6944 


6943 


6933 


6946 


6956 


2 


1 


6286 


6315 i 


6295 


6317 


6348 


3 


2 


6550 


657 0 i 


6550 


6500 


6 598 


4 


3 


692 4 


6925 1 


6910 


6927 


6939 


5 


4 


6692 


670 8 1 


6691 i 

! 


6713 


6730 


6 


c; 


6397 


6419 i 6415 ! 


6420 


6434 


7 


6 


6390 


6382 ' 6407 


6431 


64o4 


8 


7 


675 4 


6742 1 

673 3 


6722 


6727 


6750 


9 


8 


6715 


6744 


6~76l ~ ~ 


r 6762 


GAS 

MASS 

99.1 

LBM/S 


10 


U 


4484 1 


4489 4488 


4500 I 

1 


4510 


1 1 


1 


4123 1 


I 4123 4153 


4 155 


4133 ! 


1 2 


2 


42^6 ! 


i 4277 4292 


4 2 96 j 


1 4323 j 


1 3 


3 


44/5 i 4477 ! 4479 


4489 


4501 i 


14 


4 


4356 1 4353 ' 4368 


4375 


4396 


15 


5 


4161 j 4187 i 4203 


4222 


4227 


16 


6 1 


4165 i 419^. 1 4187 


4220 


4 24 5 j 


17 


7 


4365 1 439 5 | 4337 


4392 


4399 


1 8 


8 


4379 1 4373 1 4379 


4406 


4417 


GAS 

MASS 

77.4 

L8M/S 


19 


U 


i 2943 1 2957 2957 


2967 


2974 j 


20 




2756 1 2779 2779 


2601 


1 2609 j 


21 


2 


2840 f 2857 | 


2353 


2872 


i 2632 1 


22 


•a 


2943 i 2953 2953 


Z9bZ 


2970 


23 


4 I 


1 2882 ■ 2897 | 2898 


2910 


2919 


24 


5 i 


2770 2793 1 2812 


2833 


2837 


25 


6 


2773 j 2795 2611 


2821 


j 2 8A2 


26 


7 


2891 1 2901 j 


2914 


2913 


1 2921 


27 


8 


2 89 0 i 2906 


2901 


2913 


2933 


GAS 

MASS 

67.7 

LBM/S 


28 


U 


2515 j 2522 


2522 


2530 


i 2 53 6 


29 


1 


^ 2365 1 2389 


23 93 


2411 


i 2417 1 


30 


2 


2434 1 2451 


1 2451 


2462 


2471 


31 


3 


2 512 


2519 


1 2519 


2527 


2533 


32 


4 


2A6 5 


2479 


I 2480 


2490 


2493 


33 


5 


2377 


240 5 


i 2418 


2435 


2437 


34 


6 


2333 


2406 


i 2416 


2425 


2442 


35 


7 


2476 


2^81 


2492 


2496 


2499 


36 


8 


246 9 1 24 8 7 


2482 


2496 


2507 
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TABLE IV 



FLUID MASS FLOW RATE FOR VARIOUS GAS FLOW RATES, GAS FLOW 
DISTRIBUTIONS AND STEAM OUTLET TEMPERATURES 





RUN 

NO. 


DISTR. 

NO. 


STEAM TEMPERATURE 

650 625 600 5T5 550 




1 


U 


12.33 


12.54 


12.75 


12.97 


13.20 


2 


1 


11.17 


11.41 


11.58 


11.79 


12 .04 


GAS 
MASS 
X 21 . 8 
LBM/S 


3 


2 


11.63 


11.87 


12.05 


12.28 


12 .52 


A 


3 


12.30 


12.51 


12. 71 


12.93 1 13.17 


5 


4 


1 1.89 


12.12 


12 .31 


12.53 1 12.77 


6 


5 


1 1.36 


11.59 ' 11.30 ; 


11.99 j 12.21 

1 


7 


6 


11 .3 5 


11.53 1 11.73 


12.01 1 12.27 


3 


7 


12.00 


12.18 


12.36 


12.56 I 12.31 


9 


8 


11.93 


12.16 1 12.40 1 


12.63 1 12.33 


GAS 
MASS 
99.x 
L 8«/S 


10 


U 


7.9 6 


8.11 j 8.25 


3.40 




1 8.56 


11 


1 


7.32 


7.46 j 7.64 i 


7.76 


; 7.94 


12 


2 


7.60 


7.72 i 7.39 


3.02 3.20 


13 


3 


7.95 


8.09 . 8.24 


1 8.38 1 3.54 


14 


4 


7.74 


7.87 ‘ 3.03 


3.17 3.34 


i 


5 


7.3 9 


7.56 7.74 


! 7.33 i 3.02 ! 


X6 1 


1 


7.40 


7.53 7.70 


7.33 1 8 .05 1 


n 1 


7 


7.75 


T.94 1 8.07 


3.20 8.35 1 


18 i 


1 8 


7.78 


7.90 i 8.05 


3.23; 8 .33 i 

_ __ . _ . _ . _ ^ 


GAS 
MASS 
77. A 
LBM/5 


19 1 


U 


5.24 1 


j 5.34 1 5.44 


5.54 1 5.64 


20 


1 


4 .90 


5.02 ! 5.11 


5.23 1 5.33 


2x ; 


2 


5.04 


1 5.16 j 5.2o 


5.36 1 5.47 


22 1 3 


5.2 3 


1 5-33 I 


5.43 


1 5.53 5.64 


23 i 


4 


5.12 


j 5.23 


5 .33 


5.43 5.54 


2^ i 


: s 


4.92 


i -05 


5.17 


5.29 1 5.38 




6 


4 .93 


j 5.0 5 1 


5.17 


5.27 , 5.39 


26 i 


7 


5.14 


5.24 


5.36 


5.45 5.54 


27 


8 


5.13 


1 . 5.25 


5.34 


5.45 1 5.56 


GAS 

MASS 

67.7 

L6M/S 


28 U 


4.47 


i ^-56 i 


4.64 


4.72 4.81 1 


29 1 1 


A. 20 


4.32 


4 . 40 


4.50 4.59 


30 2 


4 .32 


4.43 


4.51 


4.60 i 4.69 


31 1 3 


4.46 


4.55 


4.63 


4.72 4.81 


32 ! 


1 4 

1 : 


4.33 


4.48 


4.56 


4.65 J 


4.74 


33 


t 5 


4.22 


4.3 4 


4.45 


4.55 i 


4.62 


34 


6 


4*24 


4.35 


4 . h4 


4.53 i 


4.63 


35 


7 


4 .40 


4.48 


4.53 


4.66 


4.74 


36 


a 


4.3 9 


4.49 


4.56 


4.66 1 4.76 
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TABLE V. POWER AND FLUID FLOW RATE 

FOR VARIOUS SI2E HEAT EXCHANGERS 





DISTR . 


NUMBER OF HEAT 
20 21 


EXCHA NCER 
2 3 


PASSES 

25 


UNITS 




U 


67A8 


6820 


6944 


704 3 


HP 


GAS 


1 


60 10 


6110 


62 8 6 


6433 


HP 


MASS 














121.8 


U 


11.99 


12 .1 1 


12.3 1 


12.51 


LBN/S 


LBM/S 


1 


10. 67 


10.85 


11.17 


11.43 


LBM/S 




U 


A370 


4414 


448 4 


4538 


HP 


GAS 


1 


396A 


4020 


4123 


421 1 


HP 


MASS 














99.1 


U 


7.76 


CO 

• 


7.96 


3.06 


LBM/S 


LBM/S 


1 


7.0A 


7.14 


7.32 


7.48 


LBM/S 




u 


288A 


2 90S 


2948 


2978 


HP 


GAS 


1 


265A 


2693 


2756 


2909 


HP 


MASS 














77. A 


u 


5.12 


5.17 


5.2 4 


5.29 


L9M/S 


LBM/S 


1 


A. 71 


4 .79 


4.90 


4.99 


LBM/S 




u 


2A66 


2485 


251 5 


2537 


HP 


GAS 


1 


2281 


2312 


23 6 5 


2408 


HP 


MASS 














6S.7 


u 


A .38 


4 .4 1 


4.47 


4.51 


LBM/S 


LBM/S 


1 


4.05 


4.11 


4.2 0 


4.28 


LBM/S 
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TABLE VI. POWER AND FLUID FLOW RATE FOR VARIOUS 
STEAH OPERATING PRESSURES 





DISTR. 


OPERATING PRESSURE 
400 500 600 


UNITS 




U 


6944 


6672 


6423 


HP 


GAS 


1 


6286 


6041 


5821 


HP 


NASS 












121.8 


U 


12.3 3 


12 .06 


11.82 


LBM/S 


LBM/S 


1 


11. 17 


10.92 


10.71 


LBM/S 




U 


4484 


4245 


4027 


HP 


GAS 


1 


4123 


3904 


3697 


HP 


MASS 












99.1 


u 


7.96 


7 .67 


7.41 


LBM/S 


LBM/S 


1 


7.32 


7.05 


6.80 


LBM/ S 




u 


2948 


2749 


256 9 


HP 


GAS 


1 


2756 


2 571 


2393 


HP 


MASS 












77.4 


u 


5.24 


4.97 


4. 73 


LBM/S 


LBM/S 


1 


4.90 


4.65 


4.41 


LBM/ S 




u 


25 15 


2339 


2179 


HP 


GAS 


1 


23 65 


2208 


2055 


HP 


MASS 












S6.7 


u 


4.47 


4.23 


4.01 


LBM/ S 


LBM/S 


1 


o 

f\J 

• 


3 .99 


3.79 


LBM/S 
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LAYOUT OF MARINE GAS TURBINE WASTE-HEAT STEAM GENERATOR; ONCE-THROUGH, 
COUNTER-CROSS FLOW HEAT EXCHANGER 




GAS 



Figure 1. Waste Heat Recovery Unit Layout 
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Nujmber of Heat Exchanger Passes 




General Node (m,n) 




Figure 2. 



Heat Exchanger Nodal Arrangement 



77 




( 



-UI) 




78 



Figure 3. Typical Nodal Gas and Fluid Flow 
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Figure 4. Nodal Gas and Fluid Flow for End of Pass Condition 




Figure 5. Tube Fin Configuration 




I II 

4^0 in. 0.50 in. 



Figure 6. Waste Heat Recovery Unit Tube Layout 
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Figure 8-A. Uniform Flow Distribution 
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Figure 8-B. Flow Distribution #1 
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Figure 8-C. Flow Distribution #2 
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Figure 8-D. Flow Distribution #3 
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Figure 8-E. Flow Distribution #4 



GflS FLOH DISTRIBUTION 




in o in S) in o 

cu cu — — 



X 
i — 
X 
CL 

O') 

X 

L? 



MOIJ SUD IHlOi lN20y2d 



87 



Figure 8-F. Flow Distribution #5 
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Figure 8-G. Flow Distribution #6 
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Figure 8-H. Flow Distribution #7 
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Figure 8-1. Flow Distribution #8 
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GRS MRSS FLOW RATE (LBM/S) 
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GRS MASS FLOW RATE (LBM/S) 



POWER VS. GflS MRSS FLOW RRTE 
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POWER VS. GRS MASS PLOW RATE 
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GFIS MASS FLOW RATE (LBM/S) 
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FLUID FLON RRTE VS. GRS FLOW RRTE 
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Figure 11-F. Fluid Flow Rate vs. No, of Heat Exchanger Passes 
(m =99.1 Ibm/s) 
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Figure 11-G. Fluid Flow Rate vs. No. of Heat Exchanger Passes 
(m =77.4 Ibm/s) 
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Figure 11-H. Fluid Flow Rate vs. No. of Heat Exchanger Passes 
(m =66.7 Ibm/s) 
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POWER VS. HERT EXCHRNGER PRESSURE 
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HERT EXCHRNGER PRESSURE (PSIR) 
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FLUID TEMPERATURE VS. TUBE LENGTH 
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FLUID TEMPERRTURE VS. TUBE LENGTH 
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TUBE LENGTH (FT) 
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FLUID TEMPERRTURE VS. TUBE LENGTH 
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Figure 14-E. Fluid Temperature vs. Tube Length (Distr. 
m =121.8 Ibm/s) 
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Figure 14-F. Fluid Temperature vs. Tube Length 
(Distr.'s U, 1; m = 121.8 Ibm/s) 
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Figure 15-A. Gas Temperature Distribution (Uniform, 121.8 Ibm/s) 
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Figure 15-B. Gas Temperature Distribution (Distr. 1, m = 121.8 Ibm/s) 



APPENDIX A 



Sample runs of the computer output are provided for a 
comparison of the effects of the non-uniform gas mass 
distributions on the temperatures within the heat exchanger 
and the effect on overall performance. 

Runs 1, 2 and 6 utilized a total gas flow rate of 
121.8 Ibm/s with the uniform distribution, distribution 1 and 
distribution 5, respectively. Runs 10, 11 and 15 utilized 
the same distributions but with a total gas flow rate of 
99.1 Ibm/s . Runs 19, 20 and 24; and runs 28, 29 and 33 
also utilized the same distributions but with mass flow 
rates of 77.4 Ibm/s and 66.7 Ibm/s , respectively. 
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MATRIX OF OVERALL HEAT TRANSFER COEFFICIENTS 
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MATRIX OF HEAT TRANSFER (Q) PER NODE 
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FLUID KASS FLOW BATE = 26366.6 LBM/UR 

TOTAL 0 (Bril/llRJ BASED ON THE Q 1ATRIX = 30752819.5 

U (BTU/HR) BASED ON THE HATEH/SrEAtt HATHIX = 30740245 

u (BTU/HR) BASED ON THE GAS MATRIX = 34621544.6 

THE MINI dim PINCH POINT = 3.1 DEGREES 
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ENTER NUMBER OF HEAT EXCHAN3ER PASSES DESIRED 
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flATRTX OF HATER/STEAM TEMPERATURES 
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MATRIX OF FXHAUST GAS TEMPERATURES 
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MATRIX OP 3VERALL HEAT TRANSPER COEFFICIENTS 
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tlATBIX OF BEAT TRANSFER (Q) PER NODE 
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FLUID MASS FLOW RATE = 26605. S LBH/HR 

TOTAL Q (BrU/HR) BASED ON THE Q MATRIX = 31030556.4 

Q (BTU/HR) BASED ON THE WATER/STEAM MATRIX = 31017056 

Q (BTU/HR) BASED ON THE GAS MATRIX = 26377549.7 

THE MINIMUM PINCH POINT = J.5 DEGREES 
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ENTER NUfiaER OF HEAT EXCHANGER PASSES DESIRED 
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MATRIX or rfATER/STEAM TEMPERATURES 
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MATRIX OF EXHAUST GAS TEMPERATURES 
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MATRIX OF OVERALL HEAT TRANSFER COEFFICIENTS 
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MATHIX OF HEAT TRANSFER (C>) PER NODE 
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r'LUID HASS FLOW RATE = 1 885 3 
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ENTKR EXHAUST Gh S MASS FLOW RATE (LBM/S) FOR GAS PATHS 1 THROUGH 10 
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MAThlX OF WATER/STEAH TEflPERATURES 
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HATHIX OF 3VERALL HEAT TRANSFEU COEFFICIENTS 
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MATRIX OF HEAT TRANSFER (Q) PER NODE 
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FLUID HASS FLOi< RATE = 17625.2 LBH/IiR 
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RUN NUMBER 
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MATRIX OF MATEB/STEAM rEMPERATUBES 
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MATRIX OF EXHAUST GAS TEMPERATURES 
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MATBIX OF EXHAUST GAS TEHPEBATURES 
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MATRIX OF OVERALL HEAT TRANSFER COEFFICIENTS 
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MATRIX OF iEAT TRANSFER (Q) PER NODE 
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FLUID MASS FLOW RATE = 15119.8 LBn/llR 
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ENTER NlJtlBER OF HEAT EXCHANGER PASSES DESIRED 
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MATRIX OF MATKR/STEAM TEMPERATURES 
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IIATHIX OF HEAT TRANSFER (Q) PER NODE 
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